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We investigate the properties of the Heisenberg 5* = 1 chain with bilinear and biquadratic interactions in 
a magnetic field using the Density Matrix Renormalization Group, Bethe ansatz and field theoretical consid- 
erations. In a large region of the parameter space, we identify a magnetized ferroquadrupolar Luttinger liquid 
consisting of a quasi-condensate of bound magnon pairs. This liquid undergoes a continuous pair unbinding 
transition to a more conventional Luttinger liquid region obtained by polarizing the system above the Haldane 
gap region. This pair unbinding transition is shown to be in the Ising universality class on top of a Luttinger 
liquid, leading to an effective central charge 3/2. We also revisit the nature of the partially polarized Luttinger 
liquid around and above the Uimin-Lai-Sutherland point. Our results confirm that this is a two-component liquid 
and rule out the formation of a single-component vector chiral phase. 
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I. INTRODUCTION 

A particularly interesting aspect of quantum many body 
systems is their ability to host ordered phases which emerge 
from the spontaneous breaking of one or several of the sym- 
metries of the system. In the case of spin systems, the main fo- 
cus lies on possible long-range-order (LRO) due to the break- 
ing of the SU(2) symmetry inherent to these systems. For 
5 = 1/2 Heisenberg systems, breaking this continuous sym- 
metry with a purely local order parameter implies magnetic 
ordering. However, when the spin is 5 > 1/2, there is the 
alternative possibility that the SU(2) symmetry is broken by 
a local quadrupolar order parameter, leading to spin-nematic 
phases.''^ Recent findings on NiGa2S4, a spin-1 material on a 
triangular lattice,"* indicate the possible realization of such a 
spin-nematic phase. Theoretical investigations of the bilinear- 
biquadratic S = 1 Heisenberg model on this lattice geometry 
have shown that it is possible to stabilize spin-nematic LRO 
of ferroquadrupolar and antiferroquadrupolar type depending 
on the sign of the biquadratic interaction,''-'' and that applying 
an external magnetic field leads to a remarkably rich phase 
diagram, with in particular a 2/3 magnetization plateau above 
the antiferroquadrupolar phase.'* 

In this paper we investigate the one-dimensional version of 
the model defined by the Hamiltonian 



-HSt 



n = jJ2 [cos0 S, • S,+i + sin^? (S, • S,+i) 

i 

(1) 

where J,H>0. While a variety of materials realizing 5=1 
Heisenberg systems are known,*" the biquadratic term seems 
to be more difficult to realize in nature.^ However, recent 
progress in the realization of effective spin-Hamiltonians in 
systems of ultracold atomic gases on optical lattices'* opens 
up a promising alternative route to investigate such systems in 
experiments. 

At zero field, the properties of this model are well un- 
derstood by now.""-"* For our considerations, it is helpful to 



keep in mind the following aspects of its phase diagram. Be- 
tween the two integrable points 6'tb = — 7r/4 (the Takhtajan- 
Babujian point,''*'''' TB) and 6'uls = '^/'i (the Uimin-Lai- 
Sutherland point," ULS) the chain possesses a finite Hal- 
dane gap.^^ -^ In this phase, at the so-called AKLT point 
(named after Affleck, Kennedy, Lieb, and Tasaki, Ref. 16) 
^AKLT — arctan 1 /3, the ground state is exactly known 
to be a valence bond solid, and for 9 > Saklt the spin 
correlation functions become incommensurate. ^-'^^ Between 
the ULS point and 9 = tt/2, the system is gapless and 
shows antiferroquadrupolar spin-nematic quasi long range or- 
der (QLRO).^"* For negative values of 6, between 6 = —3/4tt 
and 9 — — 7r/4 the system is dimerized and has a finite 
gap.'^"'^'^-' For the remaining values of 9, the system is in a 
ferromagnetic state. 

In comparison, the properties in a field have been inves- 
tigated less intensively. Most of the attention has been de- 
voted to the region of positive and not too large biquadratic 
interaction,-**"^' and large regions of the phase diagram re- 
main unexplored. For instance, the fate of the spin-nematic 
phase realized at zero field in the region tt/A < 9 < 7r/2''* 
remains an open issue, and the finite field properties of the 
model for negative biquadratic interactions are largely unex- 
plored, including the finite field properties of the integrable 
TB point. 

In this paper, we complete previous studies and consider 
the full phase diagram at finite magnetic fields over the whole 
range of 9 by applying the density matrix renormalization 
group method (DMRG), complemented by a Bethe ansatz 
(BA) solution of the TB and ULS points in a magnetic field 
and field theoretical considerations. We put special empha- 
sis on the possible realization of spin-nematic QLRO in the 
presence of a field by explicitly computing the spin-nematic 
correlation functions in real space. 

The scope of this paper is fourfold. First, by considering the 
full phase diagram of the system at finite magnetic fields we 
want to make precise predictions for ongoing efforts in the re- 
alization of such models in systems of ultracold atomic gases 
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FIG. 1. (Color online) (a) DMRG results for systems with L = 60 lattice sites for the magnetization of the BLBQ S — 1 chain as a 
function of {H, 9). The various colors indicate the different phases of the system, (b) Phase diagram obtained from the DMRG results for the 
magnetization, the correlation functions defined in Eqs. (3)-(5), (10) and for the central charge c as a function of [H, 9) and (c) as a function 
of (M, 9). The black line in the magnetized dimer phase indicates a crossover line at which the exponents of C^"^ and Cq^2 are both equal to 
one, see Sec. Ill C 4 (the dashed lines are linear extrapolations to the boundary of the phase and serve as a guide to the eye). The green dashed 
line below the kink transition indicates a crossover line between two different two-channel LL phases, one of them being a spin-nematic LL. 
Note that, as discussed in Sec. Ill D 2, this line is not exactly at9 — tt/A but seems to wind around this value. 



on optical lattices and for future quantum magnetic materi- 
als which eventually may be described in terms of the BLBQ 
chain. Second, we address the possibility to realize unconven- 
tional QLRO by explicitly computing the spin-nematic corre- 
lation functions in real space. Third, we consider in detail the 
phase transitions at finite field and study their critical behav- 
ior Fourth, we address the possibility to realize vector chiral 
LRO as identified previously in frustrated 5=1/2 chains in a 
magnetic field and has been proposed for the BLBQ chain in a 
magnetic field. We will demonstrate that the magnetic field 
leads to the realization of five different LL phases, and that 
these magnetic phases are connected to each other by either 
continuous phase transitions or crossovers. Our findings are 
summarized in Fig. 1 which show our DMRG results for the 
magnetization as a function of {H, 6) and the main result of 
this paper which is the complete phase diagram of the BLBQ 
chain in a magnetic field. 

The paper is organized as follows. In Sec. II we introduce 
the observables relevant for the description of the various LL 
phases. In Sec. Ill we present the complete phase diagram 
as a function of (i/, 0) by discussing our results for the Mag- 
netization (Sec. Ill A), for the central charge (Sec. IIIB), for 
the correlation functions in the single-component LL phases 
(Sec. IIIC), and for the correlation functions in the two- 
component LL phases (Sec. HID). Concerning the single- 
component LL phases, we demonstrate in Sec. Ill C 3 that a 
ferroquadrupolar LL phase is realized, and we discuss the ex- 
tension of the magnetized Haldane phase in Sec. IIIC 5. In 
Sec. Ill D 1 we demonstrate the absence of vector chiral order 
in the two-component LL phases, and in Sec. Ill D 2 we show 
that one of them is a spin-nematic LL. In Sec. IV we discuss 
in detail the transition from the magnetized dimer phase to the 
magnetized Haldane phase, which we identify to be an Ising 
transition with central charge c = 1 + 1/2 = 3/2. This sce- 
nario is further corroborated by a field theoretical treatment in 



the vicinity of the TB point discussed in Sec. V. We summa- 
rize our findings and conclude in Sec. VI. Finally, we provide 
in Appendix A and B a more detailed discussion of the BA 
solutions of the model at the TB and the ULS point at finite 
magnetic fields, respectively. 



II. OBSERVABLES AND CORRELATION FUNCTIONS 
A. Magnetization 

The properties of the model have been identified by calcu- 
lating with DMRG a number of characteristic quantities. The 
first one is the magnetization defined by 

M=i^(5f), (2) 

i 

which has been determined as a function of 9 and applied 
magnetic field H [Fig. 1(a)]. 



B. Spin correlation functions 

The second source of information comes from the behavior 
of correlation functions characterizing magnetic, spin-nematic 
and vector-chiral (quasi-)Iong-range order. We investigate 
possible algebraic decay of these coiTelation functions and 
compare the numerical values of the exponents with each 
other, the exponent with the smallest absolute value giving 
the dominant correlation function. This is of particular inter- 
est for the characterization of the gapless LL phases at finite 
field. 
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The first type of QLRO is identified by the correlation func- 
tions of the local spins. 



with AS*^ = ±2. In particular, we compute 



C trans 



(3) 
(4) 



In ID, an algebraic decay of C^"^(i,j) indicates mag- 
netic QLRO along the field, while a power-law behavior of 
C5™^(i, j) can be interpreted as magnetic QLRO perpendicu- 
lar to the field or as a quasi-condensate of magnons. 



C. Quadrupolar correlation functions 

For systems with S > 1/2, however, QLRO can also be 
identified by considering the on-site spin-nematic correlation 
functions 



CQ{^,J)^{Q^■Q3)-{Q^) -iQl), 



(5) 



where we have introduced the local quadrupolar order param- 
eter 
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(6) 



Note that only the first entry of this vector conserves S^^^^^, 
while the other ones change this quantum number by AS^ — 
1 or AS^ = 2, respectively. This leads to three components 
of the correlation functions (5): the longitudinal component 
Cq.o considering the terms conserving S^^^^^, the transverse 
component Cg.i considering the entries of (6) with AS*^ = 
±1, and the pairing component Cq^2 considering the entries 




Cc 



'Ahj) = 2 (^^{{s^st + sts- -Hsn') 

^isiSt + StSr-A{S^Y)) 



Q 



(1) 



Q 



(1) 



(7) 



CgAhj) = 7, {{St St + SfSt) [S^S- + SJS^) + h. 



'stns, 



(9) 



At zero magnetic field, these three components are identical 
due to the SU(2) symmetry of the system. At finite field, how- 
ever, the SU(2) symmetry is reduced to U(l) and the differ- 
ent components can show different behavior and character- 
ize different types of QLRO. In addition, in the presence of 
a finite field, the operators entering C^"^ and Cq^ are al- 
lowed to mix by symmetry. These correlation functions are 
thus expected to decay with the same power law and to test 
for the same type of QLRO, namely magnetic QLRO along 
the field. Similarly, Cg™^ and Cq^i can mix and test for 
magnetic QLRO transverse to the field, or for magnon quasi- 
condensation. However, the pairing component Cq 2 has no 
magnetic partner and probes possible QLRO of non-magnetic 
type. In the following we will refer to a phase at finite mag- 
netizations in which this component decays algebraically and 
dominates as a quadrupolar or spin-nematic Luttinger liquid. 
In analogy to the interpretation of (4), QLRO in Cq,2 can 
also be viewed as the quasi-condensation of 5 = 2 bound 
pairs of magnons. If, in addition, the structure factor of this 
component of the quadrupolar correlation function is peaked 
at a wavevector g = 0, we call the phase a fermquadrupo- 
lar Luttinger liquid which we sketch in Fig. 2: on a single 
site, a finite expectation value (Qi) can be envisaged as fluc- 
tuations around an axis called director^ In a ferroquadmpo- 
lar phase, all directors align parallel to each other. Note 
that the directors are perpendicular to the field. In previous 
work, ferroquadrupolar long-range-order has been identified 
for the S — I bilinear-biquadratic Heisenberg model on the 
square^^^ and triangular lattice"* at zero magnetic field. In the 
ID case, at zero field spin-nematic phases have been identified 
numerically. '"^'^^^ They are not of ferroquadrupolar type since 
the structure factor is peaked at some finite momentum. How- 
ever, in Sec. Ill C 3 we demonstrate that such a phase is real- 
ized in ID in the BLBQ chain in a magnetic field. Note that 
beyond the quadrupolar order it is possible to realize multi- 
polar (quasi-)long-range order in spin systems by considering 
bond products of spin-operators as described in Refs. 37-40. 
However, in the following we will restrict ourselves to the in- 
vestigation of possible quadrupolar order. 



FIG. 2. (Color online) Sketch of a ferroquadrupolar ordered state in 
zero magnetic field. Note that the breaking of the spin-SU(2) symme- 
try can be described by a director located on each lattice site around 
which the spins fluctuate. The parallel alignement of these directors 
perpendicular to the 2— axis is representative for the ferroquadrupo- 
lar state. 



D. Vector chiral correlation functions 

In addition to possible QLRO associated to breaking the 
SU(2) or U(l) symmetry, true LRO can be obtained by break- 
ing the parity of the system. This is tested by the first of the 
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following two vector chiral correlation functions, 
where kj = Sj x S^+i, 



(10) 



Such vector chiral order has been, e.g., identified using the 
DMRG in fmstrated S = 1/2 and S" = 1 Heisenberg 
chains In this paper, we consider a proposal of Kolezhuk 
and Vekua^''* in which a vector chiral phase at finite magnetiza- 
tions for positive values of 9 has been suggested. We address 
this issue by directly computing the correlation functions (10). 



E. Correlation exponents 

Whenever one of the above correlation functions decays as 
a power law at long distance, the decay will be described by a 
positive exponent rj according to: 



C{t,j)^\t-j\-^. 



(11) 



The exponents will be distinguished by the same indices and 
superscripts as the corresponding correlation functions: t/g"^, 
'7s™^' VQ-O^ VQ-I^ '7Q.2- If several Fourier components of the 
correlation functions decay algebraically, the exponents will 
be distinguished by an additional index. 

As pointed out above, in a magnetic field some spin and 
quadrupolar correlation functions are coupled, so that 77^"° — 
rjQ o and 77^^"^ — rjQ i. So we end up with three a priori 
independent exponents: 77^"°, rjg™'^ and 77q,2- 



F. Central charge 

The analysis of the phase diagram is complemented by 
computing the central charge c. For systems amenable to a 
description by conformal field theory this quantity character- 
izes the phase and the universality class of phase transitions. ''^ 
Using the DMRG, it can be obtained easily by computing the 
von Neumann entanglement entropy of blocks of consecutive 
sites 



5*^ = -TiQilngi, 



(12) 



where gi is the reduced density matrix of a subsystem of size 
£. In order to circumvent the oscillations which appear in the 
case of open boundary conditions (see, Refs. 36 and 44 for the 
behavior at the TB and at the ULS points at zero field), we ob- 
tain c from systems with periodic boundary conditions (PBC), 
for which Calabrese and Cardy"*^ have derived the general ex- 
pression 



5PBC ■ 



(13) 




The numerical value of c is then obtained by computing 5*^ for 
finite systems and fitting Eq. (13) to the results. 



FIG. 3. (Color online) Block entropy for a system with L = 60 
lattice sites and PBC at ^ = n/6 below the kink transition (M = 
1/3, □) and above (M = 2/3, O)- The solid lines indicate the fit 
with Eq. (13), the labels denote the central charge obtained by this 
fit. 



III. PHASE DIAGRAM OF THE BLBQ CHAIN AT FINITE 
MAGNETIC FIELDS 



In this section, we map out the phase diagram by succes- 
sively looking at the magnetization, the central charge, and 
the correlation functions. 



A. Magnetization 

The magnetization is depicted in Fig. 1(a). Three sepa- 
rate regions of magnetization can be identified. Starting from 
the dimerized phase, the magnetization of the finite systems 
grows regularly in steps of AM ~ 2/L, similar to the behav- 
ior in certain polarized 5 = 1/2 magnetic systems. 
This is rather natural. Indeed, in this parameter range (large 
negative biquadratic interactions), the system dimerizes to 
form singlet pairs, and the quintuplet crosses the triplet when 
the bilinear interaction vanishes. In this region, the first ex- 
citation of a pair of spins to cross the singlet upon applying 
a magnetic field is a quintuplet with = 2, leading to the 
step size of 2/L. In contrast to this, in an intermediate region 
around the Heisenberg point = the magnetization grows 
in steps of AM = l/L. In the region ~7r/4 < < -tt/IO 
the results for finite systems indicate a line of transitions at 
which the magnetization steps change from AM = 2/L to 
AM = l/L upon increasing 6 and keeping H — const, or 
upon increasing H at constant 9, indicating the presence of 
a phase transition. This phase transition will be discussed in 
detail in Sec. IV. By further increasing 9, we identify a kink 
anomaly developing at 6* > 7r/9 and persisting up to = 7r/2. 
Its position on the magnetization curve grows quickly up to 
the ULS point and then seems to saturate at M « 0.6, in 
agreement with previous findings concentrating on the Hal- 
dane phase (Refs. 29 and 30) and the vicinity of the ULS 
point (Refs. 28 and 31). 
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B. Central charge 



1. Exponent of the transverse correlation functions 



The central charge has been calculated throughout the 
phase diagram. It is equal to 2 below the kink, and equal to 
1 everywhere else (except along the transition line where the 
magnetization steps change from A A/ = 2/L to AA/ = 
where it is equal to 3/2, see below). Typical results are plotted 
in Fig. 3. 

When the central charge is equal to 1, the system is a single 
component Luttinger liquid, and all correlation exponents are 
controlled by a single parameter, leading to specific relation- 
ships between the exponents of algebraic correlations. 

When the central charge is equal to 2, the system is a 2- 
component Luttinger liquid, and the exponents of the alge- 
braic correlations can be obtained from the dressed charge 
matrix (see appendix B). 

So magnetization and central charge reveal the presence 
of three main phases: two single-component Luttinger liquid 
phases with magnetization steps A A/ = 2/i and AAf = 1/L 
respectively, and a two-component Luttinger liquid phase. We 
now turn to a careful investigation of correlations inside these 
phases to identify the nature of the dominant correlations. 



C. Correlation functions in the single component Luttinger 
liquid phases 



In this and the following section we base the characteri- 
zation of the phases on the values of the exponents of Cq_2, 
(^trans ^^^j (^^"s jg sufficient to consider only these corre- 
lation functions since, as expected due to the mixing at finite 
field, we find that transverse or longitudinal correlation func- 
tions decay with the same exponent, respectively. An excep- 
tion to this is the longitudinal component of the vector chiral 
correlation function, which we find to decay with an exponent 
of approx. 2 for all values of M and 0. 

The behavior of the exponents reveals important aspects of 
the phase diagram: As seen in Figs. 4(b) and (c), the trans- 
verse correlation functions decay algebraically with an expo- 
nent whose absolute value is much smaller than ?yQ,2, as soon 
as the size of the magnetization steps changes to AA/ = 1/L, 
revealing a fundamental change in the physics despite the fact 
that the central charge on both sides of the transition is c = \. 
The critical value of 9 at which this transition takes place de- 
pends on the magnetization. This scenario is reminiscent of 
the paired superfluid phase and pair-unbinding transition to 
a superfluid of single bosons identified in Ref. 48 for systems 
of hard-core bosons with correlated hopping on the square lat- 
tice. At the present, it is unclear if in this system the transition 
is of first or second order. In the case of the BLBQ chain, 
however, we identify the transition to be a continuous one as 
we discuss in detail in Sec. IV. In the following, we address 
the various aspects concerning the correlation functions and 
the physics they reveal in more detail. 



In the discussion of the phase diagram, we mainly consider 
the exponents rjg'^^^ and 77q,2- This is possible since we iden- 
tify both, from the numerical data as well as from a field the- 
ory and Bethe ansatz approach (see Sec. V and Appendix A), 




FIG. 4. (Color online) (a) Absolute value of the algebraically de- 
caying correlation functions in the ferroquadrupolar LL phase at 
M = 2/3, S = —Q.2-K. The symbols and colors indicate C''g^^ 
[Eq. (3), red line, □], Cq,o [Eq. (7), green line, ] and Cq,2 [Eq. (9), 
blue line, Ql- and C^°"= [Eq. (10), black line, vl- I" addition, 
Cg""^ [Eq. (4), magenta line, o] - expected to decay exponentially 
- is shown, (b) Value of the exponents of C^""" [Eq. (4)], of Cq,2 
[Eq. (9)], and of &p [Eq. (3)] as a function of 6 at M = 1/3. (c) 
The same as in (b) but at A/ = 2/3. The vertical lines in (b) and 
(c) indicate the position of the phase transitions and crossovers as 
depicted in Fig. 1. Note that in the high field phase aiO = 7r/4 the 
number of down spins is zero, leading to a vanishing Cq,2. In (b) 
and (c) the value of r)^^ close to the phase transition at ^ ~ — 0.27r 
is not shown since it is strongly affected by finite size effects. 
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that the exponent rjg'^^ in the single channel LL phases is the 
inverse of the one of the corresponding transverse correlation 
function, i.e., 775"^ = 1/?7q,2 in the magnetized dimer phase 
and J]^"^ — 1/77^^"^ in the magnetized Haldane phase. Note, 
however, that the longitudinal correlation functions posses an 
additional oscillating component which makes it more diffi- 
cult to obtain the numerical value of the exponent with a high 
precision. In order to obtain accurate results, we fit the expo- 
nent r?^"^ to the Friedel oscillations in the local spin density 
(5f ) using expressions obtained by bosonization, as discussed 
in Sec. V below. We apply Eq. (35) in the magnetized Hal- 
dane phase and Eq. (41) in the magnetized dimer phase by 
performing the fit only in the bulk region around the center 
of the system. We find that in the magnetized Haldane phase 
close to the transition to the ferroquadrupolar LL, finite size 
effects become predominant due to the vicinity of the critical 
point. However, in contrast to the approach used in Ref 49, 
using these expressions it is possible to obtain accurate results 
without introducing additional phenomenological fitting pa- 
rameters. Note that in the region where c — 2, there is no 
simple relation between the exponents of the various corre- 
lation functions, as discussed in more detail in Appendix B. 
However, we find this exponent to be larger than jfg^'^^ and 
VQ-2, so that we conclude that the longitudinal correlations do 
not become dominant in the two channel LL region. 

2. Oscillatory component of the longitudinal correlation functions 

We find that in the magnetized Haldane phase the oscilla- 
tory component of the longitudinal correlations decays faster 
at larger magnetizations, and substantially faster than in the 
ferroquadrupolar LL phase. The frequency of this oscillation 
depends on the value of the magnetization, and changes upon 
crossing the phase transition. It is remarkable that the wave 
vector of these oscillations in the ferroquadrupolar LL phase 
is 7r(l — A/), while at the Heisenberg point it is 27rAf. This 
can be understood in terms of the bound pairs of magnons 
populating the lattice, leading to an effective filling which is 
only half the value of the magnetization. This aspect will be 
discussed in more detail in Sec. V in the context of a field 
theoretical treatment. 



3. Ferroquadrupolar Luttinger liquid phase 

We now focus on the region of the phase diagram in which 
the magnetization steps are of size AM = 2/L at high mag- 
netizations (the light green region in Fig. 1). In Fig. 4(a) we 
present our DMRG results for the various correlation func- 
tions as obtained for systems with L = 60 lattice sites with 
open boundary conditions at Af = 2/3 and 9 — —0.2tt, a 
case which is representative for this phase. The plot shows 
that the pairing component of Cq decays slowest, consistent 
with a spin-nematic phase. Since the structure factor asso- 
ciated to this correlation function is peaked at a wave vector 
q = 0, and the central charge is found to be c = 1, we con- 
clude that the system realizes a single channel ferroquadrupo- 




FIG. 5. (Color online) Exponents Ai,2 of the power law decay of 
C^"^ and of Cq,2, respectively, at the TB point {9 — — 7r/4) as 
functions of the magnetization per site as obtained by Bethe ansatz (s. 
Appendix A). For small magnetizations C^"^ characterized by the 
exponent Ai dominates, while for large magnetization Cq^2 charac- 
terized by the exponent A2 becomes dominant. 



lar LL phase. At the same time, the transverse spin correla- 
tion function decays exponentially, showing that the single- 
spin excitation spectrum possesses a finite gap and that there 
is no quasi-condensation of magnons. This interpretation is 
further confirmed by the analysis of the one-magnon and two- 
magnon spin gap, which is presented in more detail in Sec. IV 
below. In addition, a Bethe ansatz analysis at the integrable 
TB point confirms the numerical finding of a power law decay 
in Cq 2 at finite field, while the one-magnon sector acquires a 
gap (see Appendix A). 



4. Crossover from a ferroquadrupolar Luttinger liquid to a 
spin-density-wave Luttinger liquid 

In Fig. 5 we show the exponents of the correlation functions 
Cg and Cq.2 as a function of the magnetization as obtained 
by Bethe ansatz at the TB point 6'tb = — 7r/4. We observe 
that for magnetizations Af < 0.4 the exponent of is 
smaller than one, while for larger magnetizations Cq^2 be- 
comes dominant. This is in agreement with numerical results 
for these exponents. Furthermore, we identify numerically for 
values of 6 away from the integrable TB point the existence 
of a crossover line in the low field region which we show in 
Figs. 1(b) and (c). At this crossover, the dominant correlations 
change from spin-nematic ones at large magnetizations to spin 
density wave (SDW) correlations. This is similar to a scenario 
realized by a frustrated ferromagnetic = 1/2 Heisenberg 
chain in a magnetic field, in which a crossover line divides 
spin-multipolar LL phases into a nematic and a SDW type of 
LL, as discussed in Refs. 37, 39, and 40. 



7 
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FIG. 6. (Color online) Vector-chiral correlation functions below 
(AI — 1/3, dashed lines) and above (M = 2/3, solid lines) the 
kink transition at 6 = n /6. 
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transverse correlations: 
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1.18 


1.28 


spin-nematic correlations: 










1.61 


1.38 




Bethe ansatz 


1.68 


1.20 





TABLE I. Comparison of the values for the exponents of the various 
correlation functions at the ULS point Suls ~ 7r/4 at M = 1/3 as 
obtained by DMRG and Bethe ansatz. We display the exponents of 
the non-oscillatory part, fc = 0, and the two smallest values at finite 
wave vectors ki and k2. 



1. Absence of vector chiral order 



5. Magnetized Haldane phase 

We find that the magnetized Haldane phase extends all over 
the region depicted in light blue in Fig. 1, i.e., up to values 
of 9 = 7r/2 in the high field region above the kink transition. 
For the transverse spin correlations, the exponent can be ob- 
tained with reliable accuracy and at the Heisenberg point the 
obtained numerical value compares well with previously pub- 
lished results. In the whole phase, the numerical data 



indicate f^g"® — 



trans 
S 



Interestingly, 77Q.2 ~ 'irjs 



as 



can be seen in Figs. 4(b) and (c). Both findings are in agree- 
ment with predictions from field theory presented in Sec. V. 
The exponent of Cq at low magnetizations behaves rather 
smoothly at the transition. However, at larger magnetizations 
its magnitude increases quickly when crossing the phase tran- 
sition, leading Cg 2 to decay very fast in the high field region 
of the magnetized Haldane phase. This can be related to the 
small number of down spins in that region. The same effect is 
responsible for the complete suppression of these correlations 
at the ULS point due to the absence of down spins at this point 
which is discussed in Ref. 3 1 and in Appendix A. 

Note that 77^'^"^^ jumps at the transition from the magnetized 
Haldane phase to the ferroquadrupolar LL phase. This is due 
to the nature of the phase transition and can be explained in 
terms of a field theory treatment of the transition, presented in 
Sec. VB3. 

With this we conclude the discussion of the single compo- 
nent LL phases and turn now to the behavior of the correlation 
functions below the kink transition. 



D. Correlation functions in the two-component Luttinger 
liquid phases 

In this section, we turn to the phases realized below the kink 
transition and characterize them by identifying the dominant 
correlation functions. 



The search for a parity broken phase is motivated by find- 
ings for frustrated 5=1/2 chains, in which a kink transition 
separates two single channel LLs, one of them exhibiting vec- 
tor chiral order.'*' By analogy, it has been suggested that the 
same scenario might occur in the present bilinear-biquadratic 
5 = 1 chain. ^'"^ However, as discussed in Sec. IIIB, in the 
whole region below the kink the central charge is c = 2, sup- 
porting a scenario in which two-component LL physics with- 
out vector chiral order is realized. This is confirmed by our re- 
sults for the vector chiral correlation functions shown in Fig. 6 



< 1.2 - 




FIG. 7. (Color online) Bethe ansatz results for the exponents of the 
dominant component of the correlations at the ULS point Suls = 
7r/4 in the 2 LL region below the kink. The results shown are ob- 
tained by the Bethe ansatz calculation presented in Appendix B and 
show the smallest values of the exponents of C5 and Cq.2 as a 
function of the magnetization per site. We observe a crossover be- 
tween a regime where transverse spin correlations and quadrupolar 
correlations are dominant at M ~ 0.258. Note that the exponents of 
the other components of the correlations are larger and are discussed 
in more detail in Appendix B. 
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which do not indicate the presence of parity breaking. We 
therefore conclude that vector chiral order is not reahzed and 
that below the kink transition in the whole region two channel 
LLs are realized which we will further characterize in the next 
section. 



2. Crossover from a SDW to a spin-nematic two channel Luttinger 

liquid 

As shown in Fig. 4(b), in the region 7r/9 < 9 < 7r/4 
the transverse spin correlations are more dominant than the 
quadrupolar ones, while in the region 7r/4 < 6 < tt/2 the 
quadrupolar correlations tend to be dominant. The presence 
of such a crossover line is confirmed by a Bethe ansatz anal- 
ysis at the ULS point f?uLS = 7r/4 presented in more detail 
in Appendix B. Note that the presence of the two massless 
modes leads to oscillating components at various momenta, 
which makes it difficult to obtain the numerical values of the 
exponents of the correlation functions. However, at the ULS 
point we can compare to the Bethe ansatz results, see Tab. I. 
The numerical values are obtained for M ^ 1 /3, and a good 
agreement between the DMRG and the Bethe ansatz results is 
obtained. 

As shown in Fig. 7, the Bethe ansatz results demonstrate 
that for M < 0.258 the spin-nematic correlations are domi- 
nant, while at larger fields the transverse spin correlations are. 
A further numerical analysis of the exponents around 6 = tt/A 
indicates that the crossover line is, indeed, bent. It seems to 
exist at values 9 < tt/A for M < 0.25, while for M > 0.25 
it seems to exist at 9 > tt/A, but bending back towards tt/A 
upon further increasing M. However, since the values of the 
exponents are so close, it is very difficult to identify the ex- 
act position of this crossover line numerically and therefore 
we leave this aspect for future investigations. Due to this 
crossover line, the system seems to reflect to some extent the 
behavior at zero field, where the system is in the gapped Hal- 
dane phase for 9 < ir/A, but realizes an antiferroquadrupolar 
c = 2 critical state for 9 > . Note that this spin-nematic 
LL is not a ferroquadrupolar LL since the structure factor of 
the quadrupolar correlations is not peaked at a wave vector 
q — Q, but at wave vectors which depend both on the value of 
9 and on the strength of the magnetic field H as discussed in 
Ref 3 1 and in Appendix B . 

Summarizing all these findings, we obtain the magnetic 
phase diagram of the BLBQ S ~ \ Heisenberg chain pre- 
sented in Figs. 1(b) and (c). In the next section, we will 
discuss in detail the nature of the transition from the ferro- 
quadrupolar LL to the magnetized Haldane phase at finite 
magnetizations. 



IV. PAIR-UNBINDING TRANSITION 
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FIG. 8. (Color online) Line of phase transitions between the ferro- 
quadrupolar LL and the magnetized Haldane phase (a) as a function 
of {H, 9) and (b) as a function of (Af, 9). The data points indicate 
the position at which the step size of the magnetization changes from 
AM = 2/L to AM = 1/L for systems with L = 120 lattice sites. 



a jump in the magnetization curve, showing that a continuous 
rather than a first order transition (metamagnetic transition) 
is taking place. Our findings for the magnetization curves 
indicate that an infinitesimal magnetic field at the TB point 
leads immediately to the binding of two magnons, while for 
9 > (?TB closing of the Haldane gap leads to the condensation 
of single magnons. On the other hand, for the fully polarized 
state, decreasing the field will lead to such bound states in the 
region 9 < arctan(— l/S),^"* while for larger values of 9 the 
excitations of the fully polarized state are described by single 
magnons. We therefore conclude that the line of transitions 
is located between the TB point at zero field and — 6'aklt at 
the saturation field, as shown in more detail in Fig. 8. In the 
following we provide further support that this transition is a 
continuous one. We consider the gaps 



In this section, we focus on the nature of the phase tran- 
sition between the ferroquadrupolar LL and the magnetized 
Haldane phase. By considering systems with up to L = 480 
lattice sites, we do not find any indication for the formation of 



Ai = ^ [Eo{S, + 1) + EoiS, - 1) - 2Eq{S,)] (14) 
A2 = ^ [Eo{S, + 2) + EoiS, - 2) - 2Eo{S,)] (15) 
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which have been apphed in Ref. 37 for characterizing frus- 
trated ferromagnetic 5=1/2 chains at finite magnetizations. 
Ai is a measure for the energy of single-spin excitations, 
while A2 correspondingly characterizes two-spin excitations. 
If A2 < Ai, then the lowest lying excitations are charac- 
terized as pairs of spins. We find that, after extrapolating to 
the thermodynamic limit, A2 is always zero for all values of 
the magnetization in this parameter region. This is confirmed 
by the observed algebraic decay of the quadrupolar correla- 
tion functions all over the parameter range. However, Ai is 
finite in the upper part of the magnetization curves, but be- 
comes zero at the point where the step size of the magnetiza- 
tion curves changes, and remains zero in the magnetized Hal- 
dane phase. In Fig. 9 we show results for Ai(6') after extrap- 
olating to the thermodynamic limit at three different values of 
the magnetization M = 0.3, M = 0.6 and M 0.8. In all 
three cases, to a good approximation the gap in the vicinity of 
the transition point closes linearly and remains zero after the 
transition. This supports the scenario of a line of continuous 
phase transitions. 

Next we address the universality class of the transition by 
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FIG. 9. (Color online) (a) Ai(e') at M = 0.3, M = 0.6 and M = 
0.8 after extrapolation to the thermodynamic limit. The black lines 
are linear fits in the vicinity of the critical points, (b) Example of the 
finite size scaling of the gap as a function of 1/L at M = 0.3. The 
solid lines are fits with second order polynomials and are shown as a 
guide to the eye. 



computing the central charge c on the line of phase transitions 
and in the adjacent phases. In Fig. 10 we show our results 
for the block entropy for systems with PBC and L = 60 or 
L = 120 lattice sites keeping up to to = 2000 density matrix 
eigenstates. The value of the central charge is obtained by ap- 
plying Eq. (13) for different values of 9 at fixed magnetization 
M = 0.6. At the critical point, we obtain a value close to 
c = 3/2, which is reproduced everywhere on the critical line, 
while we find the expected c = 1 in both the ferroquadrupolar 
LL phase and the magnetized Haldane phase. 

It is remarkable that the measured effective central charge 
at the transition in finite field is the same as the one at zero 
field at the TB point. In this case the transition belongs to the 
5t/(2)fe=2 Wess-Zumino-Witten-Novikov (WZWN) class,^^ 
with a unique set of scaling dimensions (and therefore expo- 
nents of the correlation functions). 

The nature of the transition in a finite field is discussed in 
some detail using field theory methods in section V. In partic- 




120 
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FIG. 10. (Color online) (a) Block entropy as obtained for systems 
with L = 120 lattice sites with PBC at 6* ^ -O.ITtt < Oc (squares), 
e « -O.IStt ^ 9c (circles) and 6 ^ -0.147r > Oc (triangles) at 
M = 0.6. The value of the central charge as obtained by fitting 
with Eq. (13) is shown, (b) Central charge as a function of 6 for 
systems with L = 60 and L = 120 sites with PBC. The data points 
are connected by a spline-interpolation as a guide to the eye. In the 
thermodynamic limit the maximum will be a sharp peak located at 
the critical point. 
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FIG. 11. (Color online) (a) Ratio R(\i — j\) [Eq. (16)] as a function 
of e at M = 0.6. (b) Value of R{30) (thick lines, filled symbols) and 
R{50) (thin lines, empty symbols) as a function of 9 at AI = 0.3, 
M = 0.6 and A/ = 0.8 obtained for systems with L = 120 lattice 
sites and PBC. The black lines are linear fits in the vicinity of the 
critical points. 



ular it is shown there that the value c — 3/2 for the central 
charge arises from an Ising degree of freedom that becomes 
critical on top of an already present critical Luttinger liquid, 
i.e. 3/2 = 1 + 1/2. This is analogous to what has been 
found in Ref. 56 for a transition between two superfluids in a 
two-component bosonic Feshbach problem, and in Refs. 57- 
61 for the case of a fermionic 5 = 3/2 Hubbard model. In the 
attractive 5 = 3/2 Hubbard model the Ising transition sep- 
arates a quasi-condensate of pairs and a quasi-condensate of 
quartets of fermions. This is reminiscent of our findings for 
the BLBQ chain in which the phase transition is connecting 
a quasi-condensate of magnons in the magnetized Haldane 
phase with the quasi-condensate of pairs of magnons in the 
ferroquadrupolar LL phase, relating a single magnon in the 
BLBQ chain to pairs of fermions in the 5 = 3/2 Hubbard 
model. 

It would be interesting to identify the Z2 symmetry of the 
Hamiltonian which gets broken at such pair-unbinding transi- 
tions. In higher dimensions, the U(l) symmetry of the system 
can be broken down partially to Z2, and this remaining dis- 
crete symmetry can further spontaneously break at the phase 



transition (see, e.g., Ref. 62). However, in one dimension, the 
continuous U(l) symmetry cannot be broken, and identifying 
the Z2 symmetry is a more difficult task. This goes beyond 
the scope of the present paper in which we focus on present- 
ing evidence in favor of such a scenario in the BLBQ chain 
at finite fields, and we leave a further characterization of this 
aspect open for future research. 

The Ising transition encountered in the models mentioned 
above can be characterized by considering particular ratios of 
correlation functions. It is shown in section V that 



(16) 



can be used as a diagnostic of the Ising transition. More pre- 
cisely, the ratio R{\i — j\) is related to a two-point function of 
an Ising disorder field /i(a;) by 



(17) 



The magnetized Haldane phase corresponds to the disordered 
phase of the Ising model, so that {^{x)) 7^ and hence R 
is expected to tend to a finite value. On the other hand, the 
ferroquadrupolar LL phase corresponds to the ordered Ising 
phase and hence R will decay exponentially. At the transition 
itself the field theory predicts R{\i — j\) ^ l/\i — j\. 

In Fig. 1 1 (a), we show our results for R at fixed value of the 
magnetization when changing 9. In the ferroquadrupolar LL 
phase, this quantity decays exponentially, while in the magne- 
tized Haldane phase it indeed tends to a finite value. At crit- 
icality, it decays ^ 1/|« — j\- Note that in this case the value 
of one or both exponents of the coiTelation functions needs to 
jump at the transition. As we will see in Sec. V, it is the expo- 
nent of (7g '^"'^ which behaves discontinuously, in agreement 
with the results shown in Figs. 4(b) and (c). In Fig. 1 1(b), we 
show the value of R{\i — j\ = 30) and R{\i — j\ — 50) as 
a function of 9 at various values of the magnetization. Even 
though we are not considering the limit of infinite separation 
between i and j, the behavior at the critical point is linear to a 
good approximation. 

Finally, we may consider the scaling of the Ising order pa- 
rameter (/i(0)) as a function of the deviation 9 — 9c from the 
critical point. As is shown in the next section, this is related 
to i?by 



lim 



R{\t~j\,9)^{ti{0))' ^\9-9c 



(18) 



Hence, the linear behavior shown in Fig. 1 1(b) is also in agree- 
ment with an Ising transition. Additional support for this sce- 
nario is given by comparing the values of the critical points 
obtained by a linear extrapolation of i?(30) and by a linear 
fit to Ai{9). From both sets of data, at M = 0.3 we ob- 
tain 9c ~ -0.2, at M = 0.6 we obtain 9c ~ -0.15, and 
at M = 0.8 we find 9c ~ —0.125. Hence we conclude that 
starting from the TB point at zero field, on the emerging line of 
phase transitions the value of c = 3/2 is kept, but the univer- 
sality class changes from SU {2)k=2 WZWN to Luttinger liq- 
uid plus Ising. This agrees with the picture emerging from the 
field theoretical treatment of this transition which we present 
in the next section. 
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V. FIELD THEORY IN THE VICINITY OF THE ISING 
TRANSITION 

Following Tsvelik*^ we can construct a field theory descrip- 
tion of the model in the vicinity of the Takhtajan-Babujian 



point Otb 
the form 



f, _ff = 0. This results in a Hamiltonian of 



'H — — LadxLa — RadxRa 



3 

imy^RgLg 



(19) 



a=l 



where La and Ra are left and right moving Majorana fermions 
and 



RbRc 



(20) 



The mass m in (19) is proportional to the deviation — Otb 
from the Takhtajan-Babujian point. The lattice spin operators 
are expressed in terms of continuum fields as 



(21) 



where x — ja^ (cq is the lattice spacing). Here AI°- are related 
to the currents M°'{x) oc J°(x), while n°-{x) are expressed in 
terms of the Ising order and disorder operators as 

n^(2;) cx a^{x)iJ,'^{x)fJ,^{x) , 
cx ^^{x)a'^{x)fi^{x) , 
n^{x) oc fi^{x)^'^{x)a^{x) . 



(22) 



We may bosonize the Majoranas 1 and 2 using 

1 



i[LiL2 + RiR2] 
Li + iL2 
Ri + iR2 



TTa 
1 

1 



■ cos V47r(c/?fl, + (Pl) 



TTa 
1 



(23) 



Here a is a short-distance cutoff. Rewriting the Hamiltonian 
(19) in terms of the canonical Bose field ^ ~ (p^ + (fj^ and 
the dual field Q = ipi^ — ipfi results in 



Ji — Hz 

iv 

2" 



-He 

'Hint 



+ Hb + 'Hint , 

L^dxL^ — Ry,dxRz] 
' 1 



K 



{dx^Y + K{dxQy 



im R3L3 
m 



cos V 47r$ 



H 



2ig 
TTa 



cos V47r<i> L3R3. 



(24) 



Here X is a function of the applied magnetic field H as 
well as the parameter 9. In order to deduce the structure of 



the ground state phase diagram of (24) we may neglect the 
marginal term Hint- This leaves us with a decoupled theory 
of an off-critical Ising model and a sine-Gordon model 
with a chemical potential equal to the applied magnetic field 
H. The latter is exactly solvable*""*'^^ and exhibits two distict 
phases. If H is less than the critical value Hc.i the model re- 
mains gapped, while it becomes critical for H > i?c,i- On the 
basis of these observations we expect altogether four different 
phases. The interaction term Hint affects the value of the criti- 
cal field Hc.i and renormalizes the parameters in H3 and He- 
In addition to the terms written in (24) further (marginal or ir- 
relevant) interactions will be generated by integrating out high 
energy degrees of freedom in the underlying lattice model. In 
particular, the marginal interaction ^marginal = dx^RsL^ 
appears to be compatible with all symmetries and therefore 
ought to be generated. While this term is unimportant as long 
as the Ising model described by H3 remains off-critical, it 
could modify the Ising transition itself. The effects of this 
interaction term have recently been analyzed in the context 
of a band-filling transition in a two-subband quantum wire by 
Sitte et al.**" Based on a 1-loop renormalization group calcu- 
lation it was suggested that ^marginal alters the nature of the 
quantum phase transition and leads to a scaling behavior that 
differs from that of H3 + He for very large length scales. 
Due to the non-scalar nature of 'Hmaiginai these length scales 
are expected to be much larger than the system sizes used in 
our DMRG computations. In the following we therefore will 

neglect "H marginal- 



A. Weak magnetic fields H < Hc,i{6) 

For weak magnetic fields H < Hc.i{0) the model (24) is 
fully gapped. The system remains unmagnetized. There are 
two phases, which are distinguished by the sign of the mass 
term in Hs-^^'^^ For 9 — 9te > we have m > and the 
Ising model described by H3 is in its disordered phase. This 
corresponds to the Haldane phase. For 9 — 9j'b < have 
TO < and the Ising model described by ^3 is in its ordered 
phase. This corresponds to the dimerized phase. 



B. Higli fields > 

When the magnetic field exceeds the critical value Hc^i{9) 
the bosonic degrees of freedom described by He enter a gap- 
less Luttinger liquid phase. For H > Hc,i{9) we end up with 
an effective Hamiltonian of the form 



"H = — [L3dxL3 ~ R3dxR3] ~ im R3L3 



V 

'2 



[dx^f + {dxQf 



(25) 



where the parameters v, v' and m depend on the magnetic 
field H and 9. This can be seen as follows. We may re- 
move the magnetic field term in (24) by the field redefinitions 

cosine term in He is oscillating in x and for sufficiently large 
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H (compared to m) drops out of the Hamiltonian J dx H at 
low energies. Carrying out a unitary rescaling of $' and Q' 
then leads to (25) (where we have also taken into account the 
renormalization of the various parameters in a strong field). 
The lattice spin operators are expressed in terms of the new 
fields as 



- cos(/3l' + 7r(l - M)x) + . . 



(26) 



where f3 > /3{Hc.i) — V^r depends on the applied magnetic 
field iy26,50 ^ jjj^j cajjjjot be easily calculated from 
within the field theory framework. In (26) we have only writ- 
ten the contributions with the smallest scaling dimensions at 
the TB point. 

The form of (25) shows that there are two phases separated 
by an Ising phase transition, which occurs when m is tuned to 
zero. We are now in a position to describe the behaviour of 
correlation functions in these phases. 



1. Haldane Phase in afield 

The case m > corresponds to the Haldane phase in a 
(strong) magnetic field. We note that the closely related case 
9 = 0, i.e. the spin-1 Heisenberg model in a field, has been 
discussed previously by several authors, see e.g. Refs [50, 51, 
63, and 68]. The Ising model described by L3, is in its 
disordered phase so that 



(a^) = , ^ 



(27) 



More precisely, the expectation value of the Ising disorder op- 
erator scales as 



{^r") cx |m| s . 



(28) 



If we consider the Ising mass as a function of the parameter 
for fixed magnetic field H, we have very close to the critical 
point 6c{H) 



m cx 10 — Oc\ 



(29) 



Using the expression given in (26) gives an exponentially de- 
caying contribution to the zz spin correlations. The leading 
long-distance behaviour is therefore due to the terms and 
sin(2/3l> + 27rMa;)^°'^', which gives 



A ^ Bcos {2nM{i-j)) ^^^^ 



Here the oscillatory contribution is always subleading as /? > 
Y^. The dominant correlations are the transverse spin corre- 
lations 

CT\^.2) - (-l)^-^(exp (^^e(x)) exp ( - 1^0(0))) 



(31) 



Hence we expect a correlation exponent for transverse spin 
correlations 



TT 1 

2/32 2' 



(32) 



The high-field phase is an attractive Luttinger liquid with 
dominant transverse spin correlations, in agreement with 
Ref 50. The long-distance asymptotics of the quadmpolar 
correlations is 



Cqaih]) ^ (exp Oye(x)) exp ( - 



(33) 



Hence we have 



l/^trans/- ■^|4 

lim J-^ = const cx 16' - 61c I, (34) 

in agreement with the numerical results shown in Fig. 11 (a). 

In order to make contact with DMRG calculations it is use- 
ful to consider Friedel oscillations in 5| for a system with 
boundaries. At the Heisenberg point 6* = we may derive 
the Luttinger liquid description of the magnetized phase of 
the open chain using a strong coupling analysis as in Refs 50 
and 51. This results in an effective spin- 1/2 Heisenberg XXZ 
chain with equal boundary magnetic fields on both ends. Stan- 
dard bosonization methods then give 



M + A- 



3in (^irMj + nip 



N+l 



sm 



\N+lJ 



(35) 



where N is the length of the chain and 

M = M 



(36) 



= 0. 



iV + 1 

We expect the form of (35) to hold also away from 
2. Dimerized Phase in a field 



Here we have to < and the Ising model described by L3, 
i?3 is in its ordered phase. Hence we have 



(37) 



As a result the contribution to the transverse spin correlations 
due to the "leading" operator identified in (26) decay expo- 
nentially, as do other contributions we have considered. The 
asymptotics of the zz spin correlator is 



COS 



{iT{l-M){i-j)) 



= (^-Jr^^cos(7r(l-A/)(^-J)). (38) 
On the other hand, the quadmpolar correlations behave as 

CQAhj) - (*~jr'"/^' ^ i^-j)-^-' ■ (39) 
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Just above Hd we have /? w -^tt, which implies that 



A, 



Ai 



(40) 



These agree with the weak-field limit of the Bethe-ansatz anal- 
ysis above. We furthermore know from the Bethe-ansatz anal- 
ysis that with increasing field the parameter grows, which 
implies that A2 decreases and Ai increases. For sufficiently 
strong fields the quadrupolar correlations become dominant. 
For an open chain we expect Friedel oscillations of the form 



sm 



^M + C- 



(^-k{1 - M)j + TTip'^j 

/32/47r ' 



jV+1 



sm 



(41) 



where for M — m/N with even m,N we find 

2(p' - M 



M = M + 



N + 1 



(42) 



3. Ising Transition 

At the transition we have rii — and the Ising model is crit- 
ical. As a result the spin correlators acquire additional power- 
law factors and become 

c'r^^hj) - i-iy-' , (43) 

Note that this implies that R{x) [Eq. 16] decays ~ 1/x, 
while comparing expressions (31) and (43) shows that the ex- 
ponent of the transverse spin correlations jumps upon entering 
the magnetized Haldane phase by a value of 1/4 independent 
of the value of the magnetization. In contrast, the exponent of 
Cq.2 changes continuously across the transition. These find- 
ings are all in agreement with the numerical results presented 
in Figs. 4(b), (c) and 11. 

In conclusion, our combined numerical and field theoretical 
analysis supports the picture that the Ising transition identified 
in the fermionic = 3/2 attractive Hubbard model finds a 
corresponding counterpart in the pair-unbinding transition of 
the 5=1 BLBQ chain at finite magnetic fields. 



VI. SUMMARY AND CONCLUSION 

To summarize, by combining extensive DMRG calcula- 
tions, Bethe ansatz and field theoretical arguments we have 
determined the complete phase diagram of the 5 = 1 BLBQ 
Heisenberg chain in a magnetic field. At finite magnetizations, 
it consists of five phases, three single component LL phases 
and two two-component LL phases. Two of the single com- 
ponent LL phases appear when polarizing the system starting 
from the dimerized phase at negative biquadratic interactions. 



At large enough fields, the LL realized in this parameter re- 
gion is a ferroquadrupolar LL, which is connected to a SDW- 
type of LL at lower fields via a crossover line. In the whole 
region, the gap to single magnon excitations is finite, and both 
LL phases are characterized by a quasi-condensate of bound 
pairs of magnons. These two single-channel LLs of pairs of 
magnons are connected by a continuous transition to the more 
standard single component LL phase of single magnons that 
appears when polarizing the Haldane phase. We determined 
the transition to belong to the Ising universality class with a 
central charge of 3/2 due to the contribution of the adjacent 
LL phases. This transition emerges at the TB point at zero 
field, showing that the magnetic field moves the universal- 
ity class from SU(2)2 WZWN at zero field to Luttinger liq- 
uid plus Ising at finite fields. The two-component LL phases 
show up for large positive biquadratic interaction (and positive 
bilinear interaction). They are separated by a magnetization 
kink from the magnetized Haldane phase and are character- 
ized by dominant incommensurate correlations of transverse 
magnetic resp. quadrupolar type. It is remarkable that they 
reflect to a certain extend the behavior at zero field. In partic- 
ular, the spin-nematic character of the LL phase identified at 
zero field in the region 7r/4 < 9 < tt/2 survives when apply- 
ing a magnetic field. It is our hope that this rich phase diagram 
will further motivate the search for experimental realizations 
of this model both, in quantum magnetic materials as well as 
in systems of ultracold atomic gases on optical lattices. 
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Appendix A: Bethe Ansatz Analysis at tlie Takhtajan-Babujian 
Point witli a Magnetic Field 

In this Appendix we describe in more detail the Bethe 
ansatz analysis of the 5 = 1 BLBQ Heisenberg chain in a 
magnetic field at the TB point 9tb = ■ At this point the 
Hamiltonian takes the form 



H 



J 

71 



(AI) 



and is known to be solvable by Bethe ansatz'"*'^^'^^ for arbi- 
trary values of the magnetic field H. The ground state is de- 
scribed in terms of the integral equation 

/02(A) = ai(A) + 03 (A) 

dfi [2a2(A - ^) + a4(A - fi)] /32(m)i (A2) 
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where 



where 



a„(A) 



1 



2n 



(A3) 



27r n2 + A2 ■ 

The integration boundary A is fixed by the condition 

£2(^)^0, (A4) 

where the dressed energy £2 (A) is a solution of the integral 
equation 



e2(A) 



.(0) 



(A) 



dfi [2a2{\ - fi) + ai{\ - ^)] €2(^1). 



-A 



Here the "bare energy" is given by 



.(0) 



(A) = 



SttJ 



ai{X)+a3{X)]+2H. 



(A5) 
(A6) 



Ground state energy and magnetization per site are 

M = l-2/ dAp2(A), 

J-A 

dX 4°) (A) p2(A). 



(A7) 



For zero field the determination of the finite-size spectrum of 
low-lying excitations is difficult because the ground state is 
made from complex solutions of the Bethe ansatz equations 
("2-strings")7'"^^ In a finite field matters are simpler and fol- 
lowing the standard analysis^^ we can establish that the finite- 
size spectrum of low-lying excited states is given by 

"(AiV)2 



2nv 
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AE 
AP ^ 2kFd 



4Z2 
2tt 



+ {Zdy + N+ + N- , (A8) 
[N+ - N- + dAN] . (A9) 



Here AA^ and d are integers and the dressed charge Z = ^{A) 
is calculated from the integral equation 

e(A) = 1 - / dfi [2a2(A ai{X - ^)] ^(m)- (AlO) 

J-A 

The result is shown in Fig. 12, which depicts the dressed 
charge Z as a function of the magnetization M. The integer 
AA^ is related to the z-component of the spin by 



~2AA^. 



(All) 



The "Fermi momentum" kp is related to the magnetization 
per site by 



fcp = 



7r(l - M) 



(A12) 



The spectrum (A9) describes a Gaussian model, which im- 
plies that the asymptotic behaviour of correlation functions 
takes the form^^ 

(O(i,x)O^(0,0)) = J2 C{d,AN,N^)e-'^''-'^^'^ 

d.AN,N± 

x{x- ivt)-'^^^ {x + ivty^^', (A13) 



2A± = 2N^ + ± Zd 



(A14) 



Which of the amplitudes C{d, AN, N^) are non-zero de- 
pends on the operator under consideration. The smallest cor- 
relation exponents for scalar operators are obtained by the 
choices 



AA^ = 0, d = 1 
d^O, AN ^1 



Ai = 2Z 



2Z^' 



(A15) 



The z-component of the spin 5| is sensitive only to states with 
AA^ = 0, while (5'^ )^ changes the total spin by 2 and hence 
couples to states with AA^ = 1. This analysis then suggests 
that the leading long-distance behaviour of correlation func- 
tions is of the form 

c'r^hj) ~ - 3)-^' cos (^(1 - Am-j)), 

CQ.2ii,j)-ii~3)-^\ (A16) 

The exponents Ai 2 are shown as functions of the applied 
magnetic field in Fig. 5. We see that for weak fields the longi- 
tudinal spin correlations dominate, while for strong fields the 
quadrupolar correlations decay more slowly. The crossover 
occurs at a magnetization per site of Mc « 0.37. 



Appendix B: Bethe Ansatz solution for 6 = Ouls = J 

Here we describe in more detail the Bethe Ansatz analysis 
of the SU(3) Uimin-Lai-Sutherland model""'^ in a magnetic 
field, which corresponds to the point Ojjls = f • For this 
value of 9 the Hamiltonian takes the form 



H 



J 



L 



(Bl) 




FIG. 12. (Color online) Dressed charge Z as a function of the magne- 
tization per site M as obtained by Bethe ansatz. Inset: magnetization 
per site AI as a function of the applied magnetic field (in units of the 
saturation field He) as obtained by Bethe ansatz (blue straight line) 
and the DMRG (red steps). 
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and is known to be solvable by Bethe ansatz for arbitrary val- 
ues of the magnetic field H. The critical properties of the 
model have been previously analyzed in Ref. [31] and we be- 
gin by summarizing the results obtained there. 



1. Ground State Properties 

The ground state is described in terms of the coupled inte- 
gral equations 



Pa(A) 

where 



b=l 



P^a\^) [ dfiK^bi)^ - fi)p,{t,) ,(B2) 



i^ii(A) = i^22(A) = -a2(A), 
pi°)(A)--27r,5„iai(A), 



(B3) 
(B4) 
(B5) 



with a„(A) defined in Eq. (A3). By virtue of the enhanced 
SU{3) symmetry of the model (Bl) in zero field the numbers 
Ma- of a spins (cr = 1, 0, —1) are good quantum numbers. 
The z-component of total spin is one of the Cartan generators 
of SU(3) and hence Af^ remain good quantum numbers even 
in the presence of a magnetic field. By definition we have 
L = Ml + Mq + and in the ground state we have 



ni = 



n2 



L 
L 



Mq + Af_i 

L 



dApi(A) , (B6) 



L 



(B7) 



-A2 



The conditions (B7) fix the integration boundaries Ai 2 as 
functions of the densities ni 2. The magnetization per site 
is 



m=^^^^=i-j: 



1 



L 

ni - n2- 



6=1 



dX pb{X) 



(B8) 



The integration boundaries v4i 2 are determined by the applied 
magnetic field through the conditions 



etiAb) = , 



(B9) 



where the dressed energies €b{X) are solutions of the coupled 
integral equations 



ea(A) = eio)(A) + ^ / dpK.tiX - (BIO) 



Here the "bare energies" are given by 



e^°\X) ^ -2TTai{X) + h , 
4^A) = /^. 



(BID 
(B12) 



1. h = 0: as a result of the enhanced symmetry the 
low-energy physics is described by the SUi{3) WZNW 
model. The central charge is c = 2. 

2. < h < hcX- the model remains in a quantum crit- 
ical phase. Universal properties are described by a 
two-component Luttinger liquid. The central charge is 
c — 2, but the symmetry is reduced as compared to 
h — 0. When h approaches hc^i the cutoff of one of the 
Luttinger liquids goes to zero. 

3. hc,i < h < hc,2' the low energy physics is described 
by a c = 1 one-component Luttinger Hquid. 

4. /ic,2 < h: the ground state is fully polarized and all 
excitations have a gap. 

In the following we concentrate on the two-component Lut- 
tinger liquid regime < h < hd- In Fig- 13 we plot the 
magnetization per site as a function of the applied magnetic 
field and the densities ni 2 as functions of the magnetization. 
We see that at the density of S*^ = —1 spins becomes 
zero. 



0.6 
0.5 



I 0.3 



0.2 
0.1 

0, 





0.1 0.2 0.3 0.4 0.5 
M 



FIG. 13. (Color online) Densities ni_2 as functions of the magne- 
tization per site as obtained by Bethe ansatz and magnetization per 
site as a function of the applied magnetic field as obtained by Bethe 
ansatz (black straight line) and DMRG (red steps). 



As a function of magnetic field there are four distinct regimes: 



2. Low-lying excitations for < h < hc,i 

As we are dealing with a quantum critical theory there are 
gapless excitations. In a finite volume L the spectrum of low- 
lying excited states scales as L^^ and is related to the operator 
content of the underlying conformal field theory^*. The finite- 
size energies and momenta of low -lying excitations can be de- 
termined by standard methods^^"^** from the Bethe Ansatz so- 
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lution with the result^' where the conformal dimensions expressed as 

2 



9Tr 1 
E{AN, d)-E, = -Y,Va{A++A-)+o{-), 

a=l 



97r ^ 



L 

a=l 



2n{nidi + n2d2) + ttANi, 

(B13) 



Af(AN,d,7Vf) = - (^Zurfi + Z2id2 ± ^^^^^^^^ -j +iVf, 

A±(AN,d, iV±) ^ I U.d, + ± ^"^"^det?^"^' )' + 



(B14) 



Here and AA^i 2 are integer numbers, 



di = — — mod 1, 



^2 = — pr— mod 1, 



(B15) 



and ui_2 are Fermi velocities of the two types of elementary 
excitations. They are given in terms of the integral equations 
(BIO), (B2) by 



(B16) 



where e^(A) are the derivatives of the dressed energies. Fi- 
nally, Zab are the elements of the dressed charge matrix 



61(^1) 62(^2) 
61(^1) 62(^2) 



where £,ab fulfil the set of coupled integral equations 



2 .A, 

£.abW = <5ab + ^ / d^l ^Qc(/i) Kcb{^J' - A) 



(B17) 



(B18) 



In Fig. 14 we plot the elements of the dressed charge matrix 
Z as functions of the magnetization. 
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FIG. 14. (Color online) Elements of the dressed charge matrix Z 
as functions of the magnetization per site A/ as obtained by Bethe 
Ansatz. 



where the exponents A are related to the finite size energies 
by 



3. Long-distance asymptotics of correlation functions 

As the critical behaviour is described by a two-component 
Luttinger liquid the asymptotic behaviour of correlation func- 
tions can be extracted from the finite-size spectrum following 
the analysis of Frahm and Korepin for the Hubbard model^''''^". 
The asymptotic behaviour of the two-point function of a local 
operator O is given by 



d,AN,N± 

^ g-2iri(nidi+n2<i2 + 5AAr 



(B19) 



A(d, AN, N+, N") = 2A+ + 2 A]; + 2A+ + 2A^20) 

For a given operator O certain amplitudes C {d, AN, ) 
will be zero due to continuous or discrete symmetries, which 
sometimes are not entirely obvious**'. 

For later use we define a number of momenta characterizing 
the oscillatory behavior of correlation functions 

Pi = 27r(7ii - ^2) , 
P2 = 7r(l - na) , 

P3 =7r(l + ni -n2). (B21) 
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4. Longitudinal Spin Correlations 



As Sj does not change the total z-component of spin only 
intermediate states for which — L ~ Ni — N2 is the same 
as in the ground state will contribute to the correlation func- 
tion. Hence the longitudinal correlations are characterized by 
quantum numbers subject to the selection rule 



AA^i + AN2 = 0. 



(B22) 



The smallest exponents are then obtained by the choices (1) 

A7Vi,2 = 0, di = ±1, d2 = 0, = 0, 

(2) AA^i,2 0, di = 0, d2 = ±1, = 0, 

(3) ANi[2 = 0,di = -d2 = ±1, iV± = 0, 

(4) AA^i^2 = 0, di = rfa = 0, N- = 0, N+ = 1, 

(5) A7Vi^2 = 0, di = da = 0, N- = 1, Ar+ = 0. 
This leads to the following form for C^"^(«, j) 



C 



'A' 



+ C2{t-jy 



cos 



(27rni(i - j)) 
[2TTn2{i - j)) 



+ C^{^-J)-^^ cos(Pi(z~j)) 



(B23) 



where 



The magnetization dependence of A^''™" is shown in Fig. 15. 
We see that the two exponents are comparable in magnitude 
but A^''^"'' < Af 

6. Quadrupolar Correlations 

Here the operator O in (B 1 9) changes the z-component of 
total spin by ±2, so that we need to consider intermediate 
states with AiVi + AN2 ~ ±2. The smallest exponents are 
then obtained by the choices 

(1) AA^i = AN2 = ±1, di = -da = ±5, N"^ = 0, 

(2) AA^2 = 2, AiVi = di,2 = iV=^ = 0, " 

This leads to the following form for Cq 

Cq.2{i.j) ^ES- COS (P3(i - i)) 



+ E2{i-j)-^Q^- 



where the exponents are given by 



(B28) 



Q.2 



A' 



Ar 



[{Z^^-Z2lf + {Z^2-Z22f] 

^11 + zl^ 



1 + det-" Z 



2 det^ Z 



det^Z 



(B29) 



The magnetization dependence of A^2 shown in Fig. 15. 



^long 

A'°"s 
A: 



2 

long 
3 



2{Zl,+Zl2). 

2(^21 + -^22) 1 
2(Zn-Z2i)' + 2(^12-^2 



(B24) 



The magnetization dependence of A'|°2^3 is shown in Fig. 15. 



5. Transverse Spin Correlations 

In the transverse spin correlator only intermediate states 
with 



ANi + AA^2 = ±1 



(B25) 



contribute. The smallest exponents are then obtained by the 
choices 

(1) AA^i = ±1, d2 = ±\, AN2 = di = 7V± = 0. 

(2) AA^2 = ±1, di = ±i, ANi = d2 = 7V± = 0. 
This leads to the following form for C^g^^^{i,j) 



<i 1.5 - 



0.1 



0.3 
M 



0.4 



0.5 



FIG. 15. (Color online) Exponents characterizing the power-law de- 
cays of the various correlators as functions of the magnetization per 
site. 



COS (P2(j-i)) 



+ L'2(i - 3) cos(7r7ii(i - 2) 



where the exponents are given by 



A trans _ r ^2 , ^2 l ^ + ^ 

Ai -[^21 + ^22] 2det2z 



A trans 
^2 



{ZI2 + ^11] 



1 + dct^ Z 
2 det^ Z 



(B26) 



(B27) 



a. Dominant Power Law Correlations 

We are now in a position to identify the dominant power- 
law correlations. In Fig. 7 we plot the magnetization depen- 
dence of the smallest exponents. We see that at low magneti- 
zations the quadrupolar correlation dominate, while for larger 
magnetizations the transverse spin correlations are seen to de- 
cay slowest. The cross-over between these two regimes occurs 
at M w 0.258. 
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